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Abstract 

i 1 Recent models of the insurance risk process use a Levy process to generalise the 

O^i traditional Cramer-Lundberg compound Poisson model. This paper is concerned 

with the behaviour of the distributions of the overshoot and undershoots of a high 

43 level, for a Levy process which drifts to — oo and satisfies a Cramer or a convolution 

equivalent condition. We derive these asymptotics under minimal conditions in the 
Cramer case, and compare them with known results for the convolution equiva- 
lent case, drawing attention to the striking and unexpected fact that they become 

_ identical when certain parameters tend to equality. Thus, at least regarding these 

quantities, the "medium- heavy" tailed convolution equivalent model segues into the 
"light-tailed" Cramer model in a natural way. This suggests a usefully expanded 

{N) flexibility for modelling the insurance risk process. We illustrate this relationship by 

comparing the asymptotic distributions obtained for the overshoot and undershoots, 
assuming the Levy process belongs to the "GTSC" class. 
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1 Introduction 



There has recently been a great deal of interest in the insurance and related literature con- 
cerning continuous time risk models, analysed under assumptions which have been found 
to be realistic for classical discrete time (compound Poisson) models. A very tractable 
continuous time generalisation of the classical model is to use a Levy process, (X t )t>o, for 
the risk process. Results for interesting quantities in this model (overshoot and under- 
shoots) can be obtained under an analogue of the classical Cramer-Lundberg condition, 
as we do in this paper (see also Biffis and Kyprianou (2010), Schmidli (1995), and their 
references for related discussions), or under a "convolution equivalent" condition on the 
tail of the Levy measure of X (see Kliippelberg, Kyprianou and Mailer (2004), Tang and 
Wei (2010), Embrechts and Veraverbeke (1982), and their references). Our aim in the 
present paper is to review the Levy setup, derive some new results for the Cramer case, 
and, further, to draw out the remarkable connection between the Cramer and convolution 
equivalent results mentioned in the abstract. This provides a very useful and practical 
synthesis of the two different approaches. 

The Levy setup is as follows. Suppose that X = {X t : t > 0}, X = 0, is a real valued 
Levy process defined on (Q, J 7 , P), with triplet (7, a 2 , Tlx), Tlx being the Levy measure of 
X, 7 e K, a 2 > 0. Thus the characteristic function of X is given by the Levy- Khint chine 
representation, E(e *) = e**( \ where 



We will consider limiting distributions of the overshoot /undershoot of the process 
above/below a high level, and the undershoot from the previous maximum, viz, the quan- 
tities X Tu — u, u — X Tii _, and u — X Tu _, conditional on t u < 00, where X t = sup 0<s<t X s , 
and r u is the first passage time above level u > 0: 



(with t u = 00 if X t < u for all t > 0). We assume either the Cramer condition, namely, 
that Ee v ° Xl = 1 for some vq > 0, or, alternatively, Ee aXl < 1 for some a > together with 
a "convolution equivalent" condition (see Section [5]) on the tail of the Levy measure of X. 
Results under the latter condition were derived in Kliippelberg et al. (2004), Doney and 
Kyprianou (2006), Park and Mailer (2008), and Griffin and Mailer (2010), which papers 
we will draw on for the relevant comparisons. The asymptotic behaviour of the process is 
quite different in these two settings, but we will make explicit some remarkable parallels 
between the associated overshoot and undershoots. 

We assume that the process X drifts to —00, modelling the situation of an insurance 
risk process with premiums and other income producing a downward drift in X, and 
claims represented by positive jumps. Thus the process X, called the claim surplus 
process, represents the excess in claims over income. We think of an insurance company 
starting with a positive reserve u, with ruin occuring if this level is exceeded by X. We 





r u = mf{t >0:X t >u}, u > 



2 



will refer to this as the General Levy Insurance Risk Model. It is a generalisation of the 
classical Cramer- Lundberg model, which arises in the special case when the claim surplus 
process is taken to be 

N t 

X t = ^2Ui-rt, r >0, (1.2) 

i=l 

where the (?7t)i=i 2,... are i-i-d. positive random variables independent of the Poisson 
process (N t ) t > . Here r represents the rate of premium inflow and Ui the size of the ith 
claim. 

The general model allows for income other than through premium inflow and a more 
realistic claims structure. Large claims, modeled by large jumps, are offset by premiums 
collected at a roughly constant rate. These claims, spaced out by independent exponential 
times, correspond reasonably well to "disasters" , of which there have been many in recent 
history. The compensated small jumps, which may be countably infinite in number, 
correspond to ongoing minor claims, with their compensator understood as the aggregate 
of premiums required to offset the high intensity of claims. The assumption lim.t_ > . 00 X t = 
— 00 a. s. is a reflection of premiums being set to avoid almost sure ruin. 

We say that "ruin occurs" if t u < 00. Interest centres on the properties of the process 
when ruin occurs, as a kind of worse case scenario, so results in this area are stated as limit 
theorems conditional on r u < 00. The assumption X t —> —00 a.s. implies P{j u < 00) < 1, 
so ruin is not certain in this model, while P(t u < 00) — » as the initial level u — > 00, so 
ruin becomes less and less probable as the initial capital is increased, as is logical. But 
P(t u < 00) > for finite u, and it is convenient to define by elementary means a new 
probability measure given by 

PM(-)=P{- \r u <oo). (1.3) 

Our convergence results will relate to this distribution. In the Cramer case they rely on 
some fundamental methods introduced by Bertoin and Doney (1994), and, like them, we 
will not need further assumptions on the tail behaviour, or otherwise, of X. 

In what follows, Section [2] lays out the basic foundational assumptions for the model, 
and Section [3] lists the marginal distributions we will need, calculated from a "quintuple 
law" . We separate the setup and results for the Cramer case in Section [4] and for the 
convolution equivalent case in Section [5] Section [6] gives a summary of the general Cramer 
and convolution equivalent results, and in Section[7]we give a diagrammatic representation 
and comparison of the results for a particular class (the "GTSC" class) of Levy processes. 
A discussion of the significance of the results is in Section [8j Proofs are collected in 
Section [9] 

2 Fluctuation Foundations 

A Levy process has certain "fluctuation" quantities associated with it as follows. Let 
L t denote the local time of X at its maximum and (L^ 1 , H t ) t > the ascending bivariate 
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ladder process of X, which is defective under the condition lim^oo X t = — oo a.s. (which 
will always obtain in this paper). The nondefective descending ladder process is denoted 
by (L^ 1 , H t ) t >Q (taking positive values). 

The defective process (L _1 ,if) is obtained from a nondefective process (C^ 1 ^) by 
exponential killing. To be precise, has an exponential distribution with some parame- 
ter q > say, and : t < Loo) = ((£~ 1 ,'H) : t < e(q)), where e(q) is independent 
of fX -1 ,?/) and has exponential distribution with parameter q. Processes (£ _1 ,"H) and 
H) are nondefective bivariate subordinators, and the marginal processes are nonde- 
fective univariate subordinators. 

Denote the bivariate Levy measure of (X -1 ,?/) by Il^-i ^(-, •). The Laplace exponent 
K,(a,b) of which satisfies 

e -«(o,6) = E(e- aL ?- bHl ;L 0O > 1) = E(e~ aC ^- mi ; e(q) > 1) (2.1) 

for values of a, b G R for which the expectation is finite, may be written 

K (a,b) = q + d L -ia + d H b+ f [ (l - e" *" 6 *) U c -i tH (dt, dx), (2.2) 

Jt>0 Jx>0 

where d^-i > and d# > are drift constants. Denote the marginal Levy measures of 
Cr x and % (equivalently of L^ 1 and H) by IL^-i and 11%, and similarly for L^ 1 and H. 
The bivariate renewal function of (L _1 ,ff) is 

poo poo 

V(t, x) = P(L- 1 <t,H s < x)ds = / e~ q8 P{L~ x <t,H s < x)ds, for t > 0, x > 0. 
Jo Jo 

(2.3) 

Its Laplace transform is given by 

e - at - bx V(dt,dx) = [ E(e- aL7l - bHs ;L OQ > s)ds = (2.4) 
>o Js>o K{a,b) 

for values a, b for which n(a, b) > 0. The marginal renewal function for H is 

poo poo 

V H (x)= P(H s <x)ds= e- qs P(H s < x)ds = lim V(t,x). (2.5) 
Jo Jo 

Analogously the renewal function of H will be denoted Vg. The tails of Tlx are 
U x (x) = U x {(x, oo)}, U x (x) = n x {(-oo, -x)}, and U x (x) = U x (x) + U x (x), x > 0, 
with similar notations Hn( x ) an d ^h( x ) f° r ^ anc ^ H- 



I I 

Jt>0 Jx 



3 Distributional Identities 

In this section we present some general formulae for distributions of the ruin variables 
(overshoot and undershoots), in terms of the fluctuation variables. From these we can 
work out the asymptotic distributions of the ruin variables. 
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The formulae are based on the "quintuple law" of Doney and Kyprianou (2009, The- 
orem 3), which gives the joint distribution of three position and two time variables asso- 
ciated with ruin. Rather than state the quintuple law in its general form, which will not 
be used in this paper, we give a formula for the marginal trivariate distribution of the 
overshoot and undershoots which can be calculated directly from it; namely, we have, for 
u > 0, x > 0, v > 0, < y < u A v , 

P (X Tu -u<x,u- X Tu _ <v,u- X Tu ^ <y,r u < oo) 
= / / U x {(v'+y',v' + y'+x)}V S (dv')\V H (u-dy')\ + d H V^(u). 

(3.1) 

Some Levy processes possess the property of "creeping" over level u, meaning that 
X Tu = u can occur with positive probability. (This cannot happen in the classical Cramer- 
Lundberg model.) By a well known result of Kesten, see Bertoin (1996, Theorem VI. 19), 

P(X Tu = u) = d H VM, (3.2) 

with the understanding that d#V#(tt) = when d# = (Vh need not be differentiable 



in this case). The second term in (3.1) is a consequence of the possibility of X creeping 
over level u. 

With the help of Vigon's (2002) "equation amicale inversee", viz 

n„(aO=/ Tt x {x + v)V Q {dv), x>Q, (3.3) 

Jv>0 

the marginal distributions of the position variables can be written down as follows. 
Marginal Distribution of the Overshoot: for u > 0, x > 0, 



P (X Tu - u<x,t u <oo)= / (n«(u - y) - Tl n (u + x - y)) V H (dy) + d H V^{u). 

J[0,u] 



(3.4) 



Letting x — > oo in this gives 

P(t u <oo)= / U H (u-y)V H (dy)+d H Vl I (u). (3.5) 



'[0. 

Marginal Distribution of the Undershoot: Define 



9v{y) ■= [ nj(z + y)V H (dz), 0<y<v. (3.6) 

Jze[0,v-y] 
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Then for u > 0, v > 0, 

P(u - X Tu _ < v, r u < oo) = / g v (u-y)V H (dy) + d H V' H (u). (3.7) 

J u—(uAv)<y<u 

For later reference we note that g v (y) is decreasing y, increasing in v, and by (3.3) we 
have 

g v (y) t Un(y) as v f oo, all y > 0. (3.8) 
Recall that for any subordinator "H, 

yn«(y)-^0 asy^O, (3.9) 



thus while it is possible that g v (0) = oo, by (3.8) and (3.9) we have 

yg v (y)^0 as y 0. (3.10) 

Marginal Distribution of the Undershoot from the Previous Maximum: for 

u > 0, < y < u, 

P{u- X Tu - < y, r u < oo) = / Ti H {u - z)V H (dz) + d H VM. (3.11) 

J u—y<z<u 

Of course this undershoot lies between and u. 

By incorporating an extension of the quintuple law which can be found in Griffin and 
Mailer (2011), similar, but slightly more complicated formulae, can be written down for 
the distributions of the ruin time and the time of the last maximum of X (i.e., the time of 
minimum surplus) prior to ruin. These will not be used in this paper, but we record them 
here for completeness. Let Gt = sup{s < t : X s = X s } be the time of the last maximum 
of X up to time t. Thus G Tu - is the time of the last maximum of X prior to ruin. We 
then have the following analogue of (3.1) for the time variables: for u > 0, s > 0, t > 0, 

P (r u - G ru - < s, G Tu _ <t)= [ n £ - s „([0, 8], (y, oo))\V(t, u - dy)\ 

J0<y<u 

+d H ^V(t,u), (3.12) 



where cL/cLw denotes left derivative. The second term in (3.12) is a consequence of the 
possibility of creeping by time t; see Griffin and Mailer (2011) for details. 
From this we obtain the 

Marginal Distribution of the Time of Minimum Surplus Prior to Ruin: for 

u > 0, t > 0, 

P(G Tu -<t,r u <oo)= / Tl n (y)\V(t,u-dy)\+d H — V(t,u). (3.13) 

J0<y<u c '- u 

Marginal Distribution of the Ruin Time: for u > 0, s > 0, 

P(r u <s)= [ [ nc-\u([0,s-t],( yi oo))\V(dt,u-dy)\+d H ^-V( S ,u). (3.14) 

J0<y<u J0<t<s a - u 
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4 The Cramer Setup and Results 



Our assumption throughout this section is the Cramer condition, namely, that 

Ee uoXl = 1, for some v > 0. (4.1) 

That this implies the process drifts to — oo, i.e., lim^oo X t = — oo a.s., and some further 
properties are summarised in the next proposition. Its proof is well known and omitted. 



Proposition 4.1 Assume (4.1). Then EX\ is well defined, and EX\ G [— oo,0), and so 
lim^ooXt = — oo a.s. Further, Ee uXl is finite and nonzero for all v G [0, uq], and Ee uXl 
increases in v for v > u ; in particular, 1 < Ee uXl < oo for v > i/q. Further, E(X 1 e uXl ) G 
(0, oo) for all v in a left neighbourhood (z/*,z/ ) ; v* < u Q , of u , and E(Xie u ° Xl ) > 
(possibly, ^(Xie^ 1 ) = +00^. 

Additionally, use of the Wiener- Hopf factorisation as in the proof of Prop. 5.1 in 
Kluppelberg et al. (2004), shows that 

E( e VoXl ) = 1 implies E{e v ° Hl \L OQ > 1) = 1. (4.2) 

Consequently, under ( |4.1 ), it follows from ( |2.2 ) with a = and b = — vq that 



q = u d H + [ (e Voh -l)U n (dh). (4.3) 

Jh>0 

Bertoin and Doney (1994) show that Cramer's estimate for ruin continues to hold for 



general Levy processes under (4.1). To state their result recall that X is a compound 
Poisson process if Ux is finite and X is the sum of its jumps. Their main result (see also 
Theorem 7.6 of Kyprianou (2005) and Section XIII. 5 of Asmussen (2002)) may then be 
stated as follows: Suppose the support of Tlx is non-lattice in the case that X is compound 
Poisson. Then 

Q 



where 



lim e U0U P(r u < 00) = (4.4) 
m* :=d H + [ he Uoh U H {dh) G (0, 00]. (4.5) 

Jh>0 

A similar result holds in the lattice case if the limit is taken through points in the lattice 



span, but to avoid repetition we will assume hence forth the non-lattice version of (4.4). 

A further useful result in Bertoin and Doney (1994) is the following version of the 
renewal theorem for Levy processes. Let 

U*{x) = f e V0V V H {dy). (4.6) 

Jy<x 



7 



Then Bertoin and Doney (1994) show that U* is a renewal function satisfying 

lim (U*(u + x)- U*(u)) = —, i£l, 

u— >oo 777,* 



(4.7) 



where the righthand side of (4.7) is interpreted as zero if m* = oo. Using (4.3) and (2.4) 
with a = 0, it is easily seen that in fact U* is the renewal function of the nondefective 
subordinator H* which has drift du* and Levy measure given by 



d H * = d H , and U H *{dy) = e Uoy U H (dy). 



(4.f 



In light of (|4.4|) and (|4.7|), it is important to know when m* < oo. From Theorem 8 

(4.9) 



of Doney and Mailer (2002b), for example, it follows that 

m* < oo if and only if EX x e v ° Xx < oo. 
Thus in the Cramer set up we will also assume 

EX ie uoXl < oo. 



We note from (4.5) that one immediate consequence of (4.10) is 

lim ue U0U Tl H (u) = 0. 

u— >co 

The following theorem is proved in Section |9j 



Theorem 4.1 Assume (4.1) and (4.10). Then we have, for x > 0, 



lim P {u) (X Tu - u<x) = l- - [ 
u^oo q J {{ 

lim P (u) (u - X Tu _ < x) 



(0,oo) 



^pdg 1 



and 



u — ^oo 

In particular, setting x = 0, 



lim ?W (u-X T ^<x) = ^ + - 

q qJo 



-l)n«(ar + dy), 
U n (y)d(e^). 



lim (X Tu 



u 



lim PM (X Tu _ 



lim P<«> (X ru _ = u) 



z/ d 



(4.10) 



(4.11) 



(4.12) 
(4.13) 
(4.14) 

(4.15) 



S 



5 The Convolution Equivalent Setup and Results 



The subexponential, or, more generally, the convolution equivalent distributions, are used 
in modeling heavy tailed data such as often occur in insurance applications; we refer to 
Embrechts and Goldie (1982), Embrechts et al. (1979), and their references, for discussion 
and properties. An important class of distributions which are convolution equivalent for 
certain values of the parameters are the generalized inverse Gaussian distributions; see 
Kluppelberg (1989). 

We briefly recap the main ideas. For background and rationale on the convolution 
equivalent assumptions see Kluppelberg et al. (2004) and Griffin and Mailer (2010). We 
say that a distribution F on [0, oo) with tail F belongs to the class & a \ a > 0, if 

lim F<y — ^ = e ax , for xG (-00,00). (5.1) 
u^oo F(u) 

Tail functions F such that F(\ogu) is regularly varying with index —a, a > 0, as u — > 00, 
are in £( a \ With * denoting convolution, a distribution F is said to be convolution 
equivalent, i.e., in the class S^ a \ a > 0, if F G & a \ and, in addition, 

F 2 *(x) 

lim := 25 Q (F) < 00, (5.2) 

where 

5 a (F) := [ e ax F(dx) 

J(0,oo) 

is the moment generating function of F. The parameter a is referred to as the index of 
the class. When a = 0, S := S^ ' is the class of subexponential distributions. 
From now on we keep a > 0. 

Distributions in the class S^ a \ for a > 0, satisfy 5 a (F) < 00 and S a+£ (F) = 00 for 
every e > 0. They are often referred to as "near to exponential" in the sense that their 
tails are "close" to the tail of the exponential distribution with parameter a. For example 
if F(x) ~ x~ 1 ~ p e~ ax for some p > and a > 0, then F G «S^"'; see Chover, Ney and 
Wainger (1973). However, note that the exponential distribution itself is not in any S^ a \ 

We can take the tail of any Levy measure, assumed nonzero on some interval contained 
in (0, 00), to be the tail of a distribution function on [0, 00), after renormalisation. With 
this convention, we say then that the measure is in or if this is true for the 
measure with the corresponding (renormalised) tail. The results of Kluppelberg et al. 
(2004) are restricted to the case where the distribution or measure is not concentrated on 
a lattice in R, which we will assume here also, with the understanding that the alternative 
can be handled by obvious modifications. Beyond this, the main assumptions are that for 
some a > 0, 

TTj(0) > 0, TT + x G S {a) and Ee aXl < 1. (5.3) 

The case when Ux is concentrated on (— 00, 0), i.e., X is spectrally negative, is rather 
easy to handle, so it is excluded throughout. (See Remark 4.6, p. 1780, of [IS].) 
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The assumption (5.3) means that we analyse the non-Cramer case. This is because if 
Ee aXl < 1 and Ee UoXT ^= 1, then < a < uq from Proposition |4.l[ But membership in 
the class means that Ee^ a+£ ^ Xl = oo for every e > 0, contradicting Ee u ° Xl = 1. Thus 
the Cramer case is mutually exclusive to the convolution equivalent case as we consider 
it here. The condition Ee aXl < 1 in (5.3) also ensures that lim^ooXt = — oo a.s. 

It's convenient to define constants 0i, /3 2 , by 



\nEe 



— ia) and (3 2 



«(0, 



-a 



Q 



(5.4) 



where recall k is given by (2.2). Under (5.3), f3\ > and < (3 2 < 1. The first is obvious 
and the second follows from Proposition 5.1 of Kliippelberg et al. (2004). 
From Kliippelberg et al. (2004, Theorem 4.2), we can deduce: 



Theorem 5.1 Assume a > and (5.3) holds. Then for x > 

lim P {u) (X Tu -u<x) = l- (3 2 e~ ax - - I (e ay - 1) Il n (x + dy) 

u ^°° q J (0,00) 



(5.5) 



Equation (5.5) compares closely with the result in (4.12), and reduces to it if we 
formally set a = uq and Ee aXl = 1, since in that case 02 = by (4.2). For x = 0, using 
(l2T2l) and ((511), rtBlij) gives 



lim P (u) (X Tu 



ad 



u 



H 



Q 



(5.6) 



as the asymptotic (conditional) probability of creeping in the convolution equivalent case. 
Compare this with the analogous result in the Cramer case given by (4.15). If d# > 0, 



then (3.2) and (5.6) imply 



lim 

w-*oo P{T U < OO) 



a 

q' 



(5.7) 



with the analogous result, with a replaced by vq, following from (4.15) in the Cramer 
case. 

The limiting distributions of the undershoots in the convolution equivalent case are 
given in Table 1 below. Their derivations using the quintuple law are explained in Doney 
and Kyprianou (2006) and carried out in Park and Mailer (2008). An alternative ap- 
proach, yielding more general results, is given in Griffin and Mailer (2010). 



6 Summary of General Results 



The general results are summarised in Table 1. In it, "Cramer Case" refers to results 
proved under Assumptions (4.1) and (4.10), while the "Convolution Equivalent Case" 



refers to results proved under Assumption (5.3), with a > 0. All the asymptotic distribu- 
tions are valid for x > 0. Recall (5.4), and note that P2 — in the Cramer case. To the 



entries in the table we can also add: 



10 



Cramer Case: 



qe 



-l/ U 



P(t u < oo) - 



as u — > oo. 



Convolution Equivalent Case: 



P(r u < oo) ~ — — Ux(u) ~ -4 n «( u )' as m oo. 

P1P2 9P2 



(6.1) 



(6.2) 



The result in the convolution equivalent case follows from Theorem 4.1 and Proposition 
5.3 in Kliippelberg et al. (2004). Note that in the Cramer case, LL^(-u) and H.%(u) are 



both of smaller order than e u ° u as u — > 00, by (4.10) and (4.11). 





Limiting distribution 
Valid for all x > 


Case I: 
Cramer case 


Case II: 
Convolution Equivalent Case 


I 


lim P (u) (X Tu -u<x) 

u— >oo 


1- 1 / (e"°y-l)Iln(x + dy) 
1 J {0,00) 


1 - p2C- aX - 1 / 

q 7(o, 


(e ay - l)U n (x + dy) 

OO) 


II 


lim P {u) (u- X Tu - < x) 

u—^oo 


is d H 1 

q q J 


f g x (y)d(e v °y) 



ad H 1 

q q J 


f g x (y)d(e a y) 



III 


lim P {u) (u~ X Tu _ < x) 

u—^oo 


v d H 1 / 

q q J 


r n n ( y )d(e v °y) 


ad H 1 / 

q q J 


f n H (y)d(e a y) 

3 


IV 


lim P {u) (X Tu = u) 

u— >oo 


v d H 

q 


ad H 

q 



Table 1: Limiting Distributions of Overshoot and Undershoots in the Cramer and Con- 
volution Equivalent Cases 



Remark: Table 1 corrects a couple of expressions in Park and Mailer (2008). The proof of 
their Theorem 3.2 is correct but there are errors in their Eq. (3.5); the simpler expression 
in Row II Column 3 of Table 1 should be used instead. Also, Eq. (3.10) of Park and 
Mailer (2008) should be corrected to that shown in Row III Column 3 of Table 1 by adding 
in the term from creeping. 

It is important to note that in the convolution equivalent case, the limiting distribu- 



tions of the undershoots are improper. Letting x — > 00 and using (3.8), we see that they 
both have total mass given by 



ad H If — ad H 1 

q — k(0, —a) 



{e ay -l)U H {dy) 



1 - ft < 1. 
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The remaining mass escapes to oo because with positive probability, as u — > oo, X can 
pass over level u as the result of an exceedingly large jump the first time it leaves a 
neighbourhood of the origin. See Doney and Kyprianou (2006) and Griffin and Mailer 
(2010) for a more detailed discussion. 

7 Example: the GTSC-Class 

In this section we use the formulae listed in Table 1 to examine the asymptotic distribu- 
tions of the overshoot and undershoots, for a particular class, the GTSC ("Gaussian- 
Tempered-Stable-Convolution") class, of Levy processes, introduced by Hubalek and 
Kyprianou (2010). For these, the marginal distributions of X\ are of generalised inverse 
Gaussian type. They generate a class of Levy processes (X t ) t > which are spectrally pos- 
itive with corresponding upward ladder height processes, (H t ) t > , being killed tempered 
stable subordinators. The downward ladder height processes are trivial subordinators 
normalised to unit drift. These processes admit some quite explicit expressions for the 
Levy measures of X and of H, and so are very suitable for calculation purposes. Our 
main aim is to illustrate similarities and differences between asymptotic distributions of 
the overshoot and undershoots for the Cramer and convolution equivalent cases. 

The GTSC class depends on five parameters g > 0, d#>0,c>0,a>0 and 
p E [—1,1). When p < we require a > 0. The Laplace exponent of X takes the 
following forms: 

$ x {0) = E[e dXl ] = -qO + d H 9 2 - cQT{-p) (a p -(a- 9f) , for 6 < a, p^0 (7.1) 

and 

*Px(0) = -q0 + d H 9 2 + c9 log -, for 9 < a, p = 0. (7.2) 

a — 9 

It can be checked that the Levy measure of X is given by 

n -(dy) = c (J^ + ^) e--dy, y>0, (7.3) 

hence LT^ e if a > and p > 0, by the result of Chover et al. (1973) referred to 
earlier. We use a for the fourth parameter of the GTSC class for this reason, so as to 
connect with the results in the previous sections. In terms of modelling, the parameters 
c, a and p determine the claim size distribution, q the premium rate and d# the noise or 
Brownian perturbation. 

The ascending ladder height process H is a subordinator killed at rate q, with drift 
d# and Levy measure 

U n (dy) = cy^ l e~ a ydy, y > 0. (7.4) 
We have H% G if, again, a > and p > 0. The Laplace exponent of H is given by 

rj) H {p) = -q + d H 9 - cT(-p) (a p -(a- 9) p ) , for 9 < a, p^0 (7.5) 
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and 



-q + d H + clog 



a 



a 



e 



, for 9 < a, p = 0. 



(7.6) 



For p G [0, 1), H has infinite activity (11% (R) = oo). For p = it is simply a (killed) 
Gamma subordinator (with drift). If —1 < p < it is a (killed) compound Poisson 
process (with drift), with intensity parameter capT(-p) and Gamma distributed jumps. 
For details of all these properties see Hubalek and Kyprianou (2010). 



From (7.1) and (7.2) we see that 

E[X 1 ]=i;' x (0) = -q<0, 



so X t — » — oo a.s. as t — > oo. Recall that the Cramer case is characterised by (4.1) 



and (4.10), while the convolution equivalent case is characterised by (5.3). For particular 



combinations of the GTSC parameters, the Cramer case is also included in the GTSC 
class. We can classify X as follows. 



Theorem 7.1 For a > 0, and X from the GTSC-class parameterised as in (7.1 ) or (7.2 ) ; 
we have the following: 

(i) if either p G [—1,0], or p G (0, 1) and d#a — q — ca p T(— p) > 0, X is in the Cramer 
case with parameter h> , where u is the unique root ofipx — on the interval (0, a); 
thus, Vq satisfies 



= ip x (vo)/vo = d H v -q- cT(-p) (a p - (a - u ) p ) , < u < a; 



(7.7) 



(ii) if p G (0, 1) and d#a — q — ca p T(— p) < 0, X is in the convolution equivalent case 
with parameter a; 



(Hi) if p G (0, 1) and d#a — q — ca p T(—p) = 0, X is in neither case. 



Proof of Theorem 7.1 (i) Recall that 

ip x is convex, tp x (0) = and ip' x (0) < 0. (7.8) 

Hence, if ipx(ot— ) > it follows that ipx has a unique zero u on (0, a) and for e > 
small enough also ipxi^o + s) < oo, giving (4.1) and (4.10). It is straightforward to check 
that the conditions in (i) indeed yield ipxi®— ) > 0. 

(ii) Next, assume p G (0, 1) and d#a — q — ca p T(—p) < 0. Then IL^ G since 
p > 0, and further 



log£[e 



ipx(ci—) = a (d#oj — q — ca p T(—p)) < 0. 



Thus Q holds. 

(iii) Finally, suppose p G (0, 1) and d#a — q — ca p T(— p) = 0. We are not in the convo- 
lution equivalent case since by the same argument as in (ii) we now have \ogE[e aXl ] = 0. 
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Nor are we in the Cramer case, since the only candidate for vq is a, by (7.8). But then it 
follows immediately from (4.5) that m* = oo and hence (4.10) fails by (4.9). □ 



It will suffice for our demonstrations to keep a > and p G (0, 1) in what follows. 
Thus we restrict to the parameter set q > 0, d# > 0, c > 0, a > and p G (0, 1). The 
parameter (3 2 defined by (5.4) here takes the form: 

/t(0, — a) —ip H (a) q — ad H + ca p T(—p) 



02 



q 



q 



(i 



(7.9) 



Thus by Theorem 7.1[ the Cramer case arises when (3 2 < and the convolution equivalent 
case when (3 2 > 0. The boundary between these regions is the surface defined by 



B ■= {p = (q, dff , c, a, p) : (3 2 = 0}. 



(7.10) 



We will denote the regions in parameter space where (3 2 < and (3 2 > by B_ and B + 
respectively. 

We will consider the asymptotics relevant to the GTSC class, using phrases like "ap- 
proach the boundary through the Cramer class (convolution equivalent class)" and vari- 
ations thereof to mean that we continuously move through points p G B_ {p G B + ) 
to a point po G B . The following observations will be used in the discussion below 
without further mention. If p° = (q°, d° H , c°, a , p°) G Bo and we approach p° through 
points p = (q,dn,c, a,p) G 5-, then one easily checks, using (7.7), that u — > a . Hence 
m* — > oo by (4.5). Of course, trivially, a — > a and (3 2 — > as p — > p° irrespective of 
whether p G B_ or p G B + . 

The following identity, which follows easily from the definition of the gamma function, 
is useful in deriving some of the formulae given in the remainder of this section; for 
< 6 < a and p G (0, 1) 



(e ey - l)y- p - L e- ay dy = -T(-p) (a p - (a - 



,-p- 1 



(7.11) 



For the probability of ruin we use the formulae in (6.1) and (6.2). In the Cramer case, 



by (4.5) and (7.4), we have 



P (r„ < oo) 



z/ (d H (a - is o y~ p - cpT(-p)) 



-UQU 



as u — y oo 



while in the convolution equivalent case, by (7.9), we have 

cq 



P (r« < oo) ~ 



a (q — adu + ca p T(— p))' 



-p-l-au 



as u — > oo. 



(7.12) 



(7.13) 



There are clear structural differences between these estimates. For example, in the Cramer 
case, the rate of decay of the ruin probability depends, through u G (0, a), on all of the 
parameters q, d#, c, a and p, while in the convolution equivalent case it depends only on 
a and p. 
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Now fix p° = (g°, d° H , c°, a , p°) G B n and approach p° through points p G £>_ or 
points p G £> + . If p G then z/ — > a , while if p G then a — > a . Thus the 



rate of exponential decay in (7.12) and (7.13) changes continuously across the boundary. 



The constant in (7.12) vanishes as we approach the boundary since m* — > oo, while 
the constant in (7.13) explodes since 02 — > 0. This, of course, is a consequence of the 
structurally different types of decay in the two cases. 

The formulae for the asymptotic distributions of the overshoot and undershoots under 
the GTSC assumption are listed in Table 2. These are calculated using the corresponding 
quantities in Table 1. 





Limiting distribution 
Valid for all x > 


Case I: 
Cramer case 


Case II: 
Convolution Equivalent Case 


I 


lim P {u) (X Tu -u<x) 


'--1 


y -p-l e - ay ^ (y-x) _ ^ dy 


1 - f3 2 e~ ax - - 

q j 


/ y-P- 1 (e- ax - e~ ay ) dy 

X 


II 


lim P {u) {u- V Tu _ < x) 

u—>oo 




- ~ I V' p ' X er av {e vay - 1) dy 

qJo 


adn c / 


i-x 

y-P- 1 (l-e- ay ) dy 


q 


+ - / 

q q J 


III 


lim P {u) {u~X Tu _ < x) 

u— >oo 


i- c 
q „ 


1 y-P- x e- av (e Voy - e u ° x ) dy 


q Jx 


y -p-l e - ay ^ ay _ dy 


IV 


lim P (u) (X Tu = u) 

u— >oo 


l - 


f -T(-p) (a"-(a-v Q y) 

q 


l - 


ca p 

q 



Table 2: Limiting Distributions of Overshoot and Undershoots for the GTSC Class 



Theorem 7.1 shows that the sign of the quantity d#a — q — ca p T(—p) distinguishes 
between the Cramer and convolution equivalent cases. In Fig.l we plot it as a function 
of a, for a particular choice of parameters (which are used for all 4 figures); namely, fix 



g = l,d H = l/2,c=l,p=l/2, 



(7.14) 



and consider the function 



f(a) := dfl-a — q — ca p T(—p) 



(7.15) 



for a G (0, 1). From Fig. 1 we see that / is negative on an interval (0,ao), the Cramer 
case, and positive on (a , 1), the convolution equivalent case, where, for the parameter 



values in fl7.14| ), a = 0.069. 

The cumulative distribution functions of the overshoot, for both cases, Cramer and 
convolution equivalent, are represented in the next diagram as a varies, for the same 
choice of parameters as in Fig 1. Recall that varying a is a means of varying the claim 
size distribution. Even though increasing a beyond an results in changing the model 
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Figure 1: The function / defined in (7.15) for values < a < 1. Parameter values are 
q=l, d H = l/2,c=l,p=l/2. 



from one in which the Cramer condition holds to one in which it doesn't, the overshoot 
distribution transitions smoothly throughout the entire range of a. A similar phenomenon 
can be observed by varying the other parameters. 



a =0.1 




Figure 2: Cumulative distribution of the overshoot as given by Row I of Table 2, in the 
Cramer case for values of a = 0.07,0.08,0.09,0.10 (full lines), and in the convolution 
equivalent case for values of a = 0.03,0.04,0.05,0.06 (dotted lines). Parameter values as 
for Fig. 1. 

The formulae for the overshoot in Row I of Table 2 are obtained by direct substitution 
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of (7.4) into Table 1. For computational purposes these formulae can be expressed in 



terms of the incomplete Gamma function, and easily calculated for specific values of the 
parameters. 

As the boundary is approached, either through B_ or B + , the formulae agree in the 
limit, for every x > 0. Figure 2 illustrates this in the case that the boundary is approached 
by varying a only. Setting x = gives the asymptotic (conditional) probability of creeping 
over the level u. The resulting formulae agree with those in Row IV of Table 1: In the 
Cramer case 



lim P (u) (X r 



u 



Q Jo 

i + -r( 
q 



p) {of 



\)y- p - x c' av dy 



[a 



(7.16) 



where the second equality holds by (7.11) and the third by (7.7), while for the convolution 



equivalent case, by the same means but using (|7.9|), we find that 
lim P {u) (X Tu 



u 



1 



ccr 



02 + — T(-p) 
1 



dncn 



(7.17) 



Comparing (7.16) and (7.17), note that in both cases there is creeping in the limit 



if and only if d#> 0, which is in turn equivalent to X having a Gaussian component 



(namely, from (7.1) we see that ax = v / 2dff)- One obvious difference between the two 



cases is that in the Cramer case changing any parameter changes vq and hence (7.16) 



while (7.17) does not depend on p and c (we saw a similar effect for the ruin probability). 
Next, for the behaviour when x — > oo, we can calculate, in the Cramer case 



1 - lim P {u) (X Tu -u<x)~-(( 



,-i 



a 



-i 



x - P -i e - ax 



x — > oo, 



while in the convolution equivalent case 

1 - lim P (u) (X Tu - u<x) 



foe 



as x — y oo. 



The structural difference between the two cases is largely analogous to that seen for 
the ruin probabilities, but curiously with the roles reversed. In this case, approaching 
the boundary through the Cramer class causes the coefficient to explode, whereas the 
coefficient vanishes if the approach is through the convolution equivalent class. 

This analysis shows that, while the formulae in the two cases are structurally different, 
especially in their asymptotic behaviour, nevertheless they segue continuously into each 
other across the boundary Bq. This effect is apparent in Fig 2. 

Similar analyses can be given for the undershoots. 

For the undershoot itself, from (3.3), (3.6) and the fact that Vg(dy) = dy since H is 
unit drift, we obtain 

9x(dy) = -U x (y)dy = -U H (dy). 
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Figure 3: Cumulative distribution of the undershoot as given by Row II of Table 2. Values 
of a and parameter values as for Fig. 1. 



Integrating by parts in Row II of Table 1, then yields the results in Row II of Table 2: 



lim P (n) (n-X Tu _ <x) 
in the Cramer case, and 



' /,,d ' rl ' 1 - -P- I e~ ay (e Voy - 1) dy (7.U 



+ 



y 



lim PW {u~X T ^<x) = ^ + - 



Q QJo 



f y^ 1 (1 - e- av ) dy 
Jo 



(7.19) 



in the convolution equivalent case. Just as for the overshoot, there is a continuous tran- 
sition across the boundary in the formulae for every x > 0. When x = they reduce 
to 



lim P(«> (X Tu . 



u 



and lim P (u) {X Tu „ = u) 



ad 



H 



(7.20) 



respectively. When x — > oo, the right hand side of (7.18) converges to 



**** - £r(- p ) (a" - (a - z/ ) p ) = 1, 
q q 



by (7.7) and (7.11), while the right hand side of (7.19) converges to 

ad.H ca p T(—p) 



1-/32 



(7.21) 



(7.22) 



by (7.9) and (7.11 ). Since /3 2 > in the convolution equivalent case, this last expression is 



in (0, 1). Hence in the convolution equivalent case, the undershoot converges under 
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to an improper distribution, which has an atom at oo of mass This atom vanishes as 
we approach the boundary through B + . 

Just as for the overshoot, another difference between the cases is in their asymptotic 
behaviour, namely, in the Cramer case we see 



1 - lim P {u) (u - X Tu _ < x) 



while in the convolution equivalent case 

l-fa- lim P {u) (u - X Tu „ < x) ~ 



PQ 



-x 



as x — y 00. 



(7.23) 



(7.24) 



As the boundary is approached through the exponential factor in (7.23) vanishes, 
but the coefficient explodes. In contrast, the estimate in (7.24) behaves smoothly at the 
boundary. 

For the undershoot from the previous maximum, the formulae in Row III of Table 2 
follow from Row III of Table 1 after an integration by parts and substitution of (7.4). 
Just as for the overshoot and the undershoot, there is a continuous transition across the 
boundary in the formulae for every x > 0. When x = 0, upon using (7.11), they reduce 
to 



lim PW (X Tu _ = u) = and lim P« (X 7 

u— >oo u ' q u— >oo 



ad 



u 



H 



Q 



(7.25) 



respectively. Precisely the same asymptotics as x — > 00 apply as in (7.21) and (7.22). 
As for the undershoot, in the convolution equivalent case, the undershoot of the last 
maximum before t u converges under to an improper distribution, with an atom at 
00 of mass /?2- This atom vanishes as we approach the boundary through B + . Finally 
(7.24) continues to hold in the convolution equivalent case if A Tu _ is replaced by A ru _, 
while in the Cramer case 



1 - lim P (u) (u - X Tu „ < x) ~ -((a - uo)- 1 - oT^aT^V^-^, 



as x — y 00. 



8 Discussion 

The quintuple law, referenced in Section |3j is the result of a deep analysis into Levy 
process theory, but the derivations from it of distributional identities such as those in 
Section [3j and, from them, the applications to limit laws such as those in Sections [4] and 
[5j are straightforward. In fact, once given the results in Section [3j our paper is virtually 
self-contained, requiring in addition nothing more than some standard renewal theory and 
the dominated convergence theorem for the proof of Theorem |4.1| The proof of Theorem 



5.1, also based on the quintuple law (see Kluppelberg et al. (2004)), requires some extra 
techniques from the theory of convolution equivalent distributions, but these are standard 
and easily applied. 
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Figure 4: Cumulative distribution of the undershoot from the previous maximum as given 
by Row III of Table 2. Values of a and parameter values as for Fig. 1. 



Consequently, we can put forward the methods exemplified in this paper as a simple 
and clear way to set out and derive limit laws for overshoots and undershoots in both 
the convolution equivalent and Cramer formulations. On the other hand, the method is 
not a panacea for all problems. Inspection of Eq. (3.1) for the joint distribution of the 
3 positions variables, and of Eq. (3.12) for the joint distribution of the 2 time variables, 
suggests that our methods will not carry over easily to these situations. We have useful 
results available in these cases, but they are derived from a deeper path analysis of the 
Levy process and will be presented separately. 

As we illustrated in Section [7], the formulae listed in Table 1 can be quite easily 
calculated in special cases, and it would be very useful to develop computational and/or 
simulation methods to deal with more general cases. 
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9 Proofs 



Proof of Theorem 4.1 Assume (4.1) and (4.10). Subtracting (3.4) from (3.5) and 
ing by P(t u < cx 

P (u) {X Tu -u>x) 



dividing by P(r u < oo) gives, for x > 0, 

1 



P(t u < oo 



q- 1 u m*e uou 



U n (u + x - y)V H (dy) 



[0,u] 



U n (u + x-y)V H (dy) (by (4.4)) 



[o,«] 



f e v °^U H (u + x-y)U*(dy) (by Q)) 
i[o,«] 



Thus, integrating by parts, 
pM(X Tu -u>x) ~ 



g 1 u m : * 



q- l v m* / e^n w (x + 2 /)d 2/ (U*(u) - U*(u - y)) . (9.1) 

'[0,u] 



q- 1 u m*e' /0U U n {u + x)U*{u) - 

q~\m* [ (U*(u) - U*(u - y)) d y (e»°m H (x + y)). (9.2) 

J[0,u] 



For later reference, note that for any x > 0, 



/ y\d y (e^U n (x + y))\ < [ 

J[0,oo) J[C 



ye Vo m H {x + dy) 



yu e^U n (x + y)dy- 
= [ [ iy ze UoZ U n (x + dy)dz+ [ ye Uoy U n (x + dy) 

J[0,oo) J[0,y) J[0,oo) 

< 2 I ye U0V Il n (x + dy) < oo, (9.3) 



by (4.5) since m* < oo. 

First assume d# > 0. Then d#* = dn > by (4.8), and so by Prop III. 1 of Bertoin 
(1996), U*(y) < ciy for all y > 0, for some constant c\ > 0. Thus the first term on the 
RHS of (9.2) tends to as u — > oo, by ( |4.11[ ). Next by subadditivity of U*, we have 



< U*{u) - U*(u — y)< U*(y) < c x y, for all < y < u. 



(9.4) 



Thus by (|9.3|) we may apply dominated convergence to the second term on the RHS of 

d y (e™TL u (x + y)) 



( [9~2| . It then follows from ([D} and Q4.11fl that 
lim P {u) {X, 



u > x) 



q 1 u m* 



q 1 "o 



y_ 

[0,oo ) vm * 

e vo m H (x + y)dy 



[0,oo) 



q- 1 [ (e^-l)U H (x + dy), 

J\0,oo) 



'[0,oo) 

21 



for all x > 0, proving (4.12) when d# > 0. 



Now assume dn = 0. In that case we can only deduce from Prop III. 1 of Bertoin 
(1996) that for each y > 0, there is a constant c 2 = c 2 (?/o) > such that U*(y) < c 2 y 



for y > y , and so (9.4) has to be modified accordingly. To account for this we need to 



consider separately the cases x > and x = 0. The first term on the RHS of (9.2) tends 



to 0, as before, in either case. When x > we can change variable in the second term to 
get 

(U*(u) - U*(u -(y- x))) d y (e^U H (y)). 



[a;,a;+u] 



By subadditivity, the integrand is bounded above by U* (y — x) < U* (y) < c 2 y for some 



c 2 = c 2 (a;), if y > x > 0. Since y\d y (e u ° y Tl-}{(y))\ is integrable by (9.3), we can apply 
dominated convergence to obtain 



lim P {u) (X Tu -u> 



x 



-q- 1 u m*e- U0X 



[a;,oo) 



y-x 



d y (e^Tl n (y)) 



q- 1 [ (e^-l)U H (x + dy), 

J[0,oc) 



just as before. It remains to deal with the case d# = and x = 0. But then (4.12) follows 



immediately from (3.2) and (4.3) 



Next we prove (4.13) and (4.14) for x > 0. We first observe that for any x > 



e^U H (y)d y (U*(u)-U*(u-y)) 

(x,u] 

= e^ u Tl H (u)U*(u) - e^Tl H (x) {U*(u) - U*((u - x)-)) 
l {x <y<u } (U*(u) -U*(u- y)) d y {e v °m H {y)). 

(0,oo) 

Now for some constant c 2 = c 2 (x), U*{u) — U* {u — y) < c 2 y for y > x > 0, as observed 



above. Thus, as u — > oo, the first term tends to by (4.11), while 2/|d y (e^ oy ri-^(2/)) | is 



integrable on (0, oo) by (9.3). So by dominated convergence 



e™n n (y)d y {U*(u)-U*(u-y)) 



(x,u] 



X 



— > — I — - ) e V0X U n (x) - 



(x,oo) 



mr 



- d y (e^U n (y)) 



mr 



e^U H {y)dy. 



(9.5) 



(x,oo) 
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Now consider (4.14) for x > 0. Subtract (3.11) from (3.5) to obtain, for < x < u, 

1 



P (tt) {u - X Tu _ > x) 



P (r„ < oo) 



~ q- 1 u m*e U0U 



[0,u-x) 



U n (u-y)V H (dy) 
Ti H (u-y)e-™U*{dy) 



[0,u-x) 

= q-'vom* [ e v °m n (y)d y (U*(u)-U*(u-y)) 

J (x,u] 

Q J(x,oo) 



by (9.5). Since 



f e^U H {y)dy = - f (e™ - 1) U H (dy) = 1 - ^£ 

9 V(0,oo) 9 </(0,oo) 9 



by (4.3), (4.14) then follows. 



To prove (4.13) for x > 0, subtract (3.7) from (3.5) to obtain 



P(u- X Tu _ > X, T u < oo) 



n w (« - y) Vh (dy ) - / (tt - y ) V H (dy) . 

[0,tt] >/ yd[u~ x, u] 

(9.6) 



The first term on the RHS of J9~6] ), divided by P{r u < oo), is just P (u) (X Tu - u > 0), by 

(9.7) 



the first equality in (9.1). Applying the already proven (4.12) with x = 0, we get 

vod H 



lim / U n (u - y)V H (dy) = lim P<«> (X Tu -u>0) = l 

t_>0 ° </[0,u] 



The second term in (9.6), when divided by P(r u < oo), is 
1 



p{t u < oo) y [u _ X)1t; 



^(«-y)e-^C/*(dy) 



Q J[o,x 



g x (y)e^d y (U*(u)-U*(u-y)). 



(9.8) 



Take 5 G (0,x) and write Jj 0x j = Jj 0iJ j + C 5a .]. For the integral over [0,5] we have 



/ g x (yy oy d y (U*(u) -U*{u- y)) < [ U n (y)e^d v (U*(u) - U*(u - y)) (9.9) 

J [0,6] J[0,S\ 



by Q. On the RHS of Q, write L s] = L } - L } , then use QsTTJ and fl4l2| |, both 



with x = 0, and (9.5) with x = 5, to see that, as u — > oo, the RHS of (9.9) converges to 
1 



(e 1 ™ - 1) U H (dy) e^U H (y)dy = —[ e^U H (y)dy. 

^m* 7 [0,oo) m J(S,oo) m J[0,S\ 



(9.10) 
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For the integral over (S,x] in (9.8), integration by parts gives 



g x (y)e™d y (U*(u)-U*(u-y)) 

g x (x)e v ° x (U*(u) - U*((u - x)-)) - g x (5y oS (U*(u) - U*((u - 5)-)) 
(U*(u)-U*(u-y))d y (g x (y)e™). 



(8,x] 



(9.11) 

We can apply bounded convergence to this integral since \d y {g x {y)e uoy ) | is a finite measure 
on (S,x\, and the integrand is bounded because U*{u) —U*(u — y) < U*(y) < U*(x) for 
< y < x. Thus 



g x (y)e™d y (U*(u)-U*(u -y)) 



(S,x] 



X 

m' 
1 

m* 



m* J(s,x] ^ m ' 

g x (y)e v ° y dy, 

(S,x] 



(9.12) 



after another integration by parts. Combining (9.6)-(9.10) and (9.12), and letting 5 — > 



proves (4.13) for x > 0. 



To complete the proof we must show (4.13) and (4.14) hold for x = 0, that is, the 



last two equalities of (4.15) hold. By (5.2) of Griffin and Mailer (2011), for any Levy 
process, P(X Tu ^ < u = X Tu ,r u < oo) = 0. Since X Tu - <«on {r u < oo}, it follows that 
P{X Tu = u,r u < oo) < P(X Tu _ = u,t u < oo). On the other hand X Tu - < X Tu _ < u on 
{t u < oo}, thus we may conclude that 

P<»> (X Tu =u)< PW (X Tu _ =u)< PW (X Tu _ = u) . 



Since (4.12) has already been proved for x = 0, we have 

vod H 



lim P (u) (X Tu =u)> 



Q 



(9.13) 



while for every x > 0, 

limsupP (u) (X Tu _ =u) < lim P {u) (u - X Tu _ < x) 

u— >oo it— >oo 

= ^ + i rn^)d(e-) 



(9.14) 



by (4.14). Letting x — » completes the proof. 



□ 
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